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HYDRODYNAMICAL VORTEX ON THE PLANE 

1. INTRODUCTION 

Hydrodynamics has always been a source of nonlinear equations integrable by the 
inverse scattering method. To mention a few, the Korteveg - de Vries, Kadomtsev - 
Petviashvili, Davey - Stewartson and Bendjamin - Ono equations all have hydrodynamical 
origin. As a rule, these equations can be obtained by direct application of procedures of 
\ the multiscale decomposition type to the equations of hydrodynamics. In other words, 

virtually all of them arise as a result of certain reductions of the initial system of nonlinear 
hydrodynamical equations, the latter having rather few known exact solutions (one should 
mention here the Landau solution for the submerged stream [1] , and a rapidly developing 
approach to the integration of nonlinear equations based on the group theory [2]). It is 
C/} \ natural then to search for completely integrable partial cases of this system, since such 

cases, if found, would allow regular procedures for constructing the solutions. This idea 
has apparently first been realized in [3], where a Lax pair for the Euler equation for ideal 
incompressible fluid was found. A more special case - the vortex equation - was studied 
in [4]. 

In this paper we consider vortical solutions of the 2D Euler equations, which have nu- 
merous applications, in particular, in geophysical hydrodynamics [5,6] and plasma physics 
[7]. The existence of a Lax pair (although, in a "weak" sense) for the problem allows to 
apply a range of tools from the theory of completely integrable systems. 
\q ', Notice that the structure of the pair under consideration is rather special - since both 

equations of the pair are of the first order in the auxiliary function, no scattering problem 
occurs. This leads, in particular, to solutions containing arbitrary functions and certain 
other peculiarities. 

The structure of the paper is as follows. In section 2 we give the necessary pre- 
liminaries from the hydrodynamics and the theory of vortices. Section 3 is devoted to 
solving a simplified vortex equation, section 4 to constructing exact 2D solutions of the 
vortex equation, both partial and more general, obtained by application of the Darboux 
transformation. 

2. PRELIMINARIES 

Let us consider the flow of an ideal incompressible two-dimensional fluid in a domain 
D G M 2 . Let v = v(x, y) = (v\, V2) be the velocity field, p be the pressure, and p be the 
density of the fluid. Then the following system of equations on the velocity field holds [1] : 

div v — 0, 

(2.1) 
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dv Vp 

m + {{JV)V = 

The first equation in (2.1) expresses the conservation of mass, the second one is the Euler 
equation. 

This system should be provided with boundary conditions. It is natural to assume 
that if the domain D is bounded, then 

(v,n) = at (x,y) e <9B, (2.2) 

where n is orthogonal to the boundary surface, which means that the fluid can not spill 
over the boundary, and that 




v — > at V^ 2 + V 2 -> oo, (2.3) 

if the domain is unbounded. 

As is well known, the system (2.1) have an infinite series of the first integrals [8]. In 
particular, the quantities ( k = 1,2 ) 

■&-«T )<to * (2A) 

are first integrals. 
Using the identity 

-Vf 2 = v A rot v + (vV)v, 
we find from the second equation in (2.1) that (dw = dp/ p) 

_ v 2 . 

— v A rot v — —\/(w H ). 

dt v 2 1 

Applying the rot operation to this equality, we obtain, 

drot v r . . _ „. 

— = rot \v A rot v \. (2.5) 

This equation is called the vortex equation. It shows, in particular, that in the stationary 
case the field of a vortex commutes with the velocity field. Assuming, further, that v\ = 
—ty y , i>2 = ^x, where \I/ = ty(x,y,t) is the stream function 1 , we find that rot v = kAty, 
where A = d 2 x + d 2 y is the Laplace operator, k is the unit vector in the direction of z 
- axis. Setting Q = A\l/ and defining the Jacobian (or hydrodynamical bracket) by the 
relation 

J{w 1 ,w 2 ) = ^ = {w u w 2 } H = w lx w 2y - wi y w 2x , (2.6) 



1 Given the velocity field, it is always possible to restore the stream function: 



^(x,y,t) = / v 2 dx - vidy, 

J (xn.Vn) 



'{xo,Vo) 

where (xq, yo) G D is a reference point, and the value of the integral does not depend on the contour of 
integration. 
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we obtain the following nonlinear equation for vorticity Q 2 : 



Q t = {Q,^} H . (2.7) 

It can be checked by direct calculation that it is possible to re-write this equation in the 
following hamiltonian form [9], 

— = {n, H} PB , n=-J (W) 2 d* dy, (2.8) 
where 7i is the Hamiltonian, and the Poisson bracket on smooth functionals is defined by 

= J il U xd y. (2.9) 

To state the problem in a closed form, it is now necessary to restore the stream function 
from vorticity. The corresponding equation is reduced to an interior Dirichlet problem for 
the Poisson equation (F is a given function), 

A* = n(x,y,t), V ldo = F {l). (2.10) 

Its solution has the form 



*(x,y,t) = A~ 1 n= [ [ G (x-x , ,y-y')Q(x , ,y , ,t)dx'dy , + [ ^-F (l)dl, (2.11) 
J Jo Jdo on 

where Go(x,y) = (l/27r)lnr + go(x,y) is the Green function for the Dirichlet problem, 
r = \/ x 2 + y 2 , and go(x,y) is a harmonic function such that Gq(x, y)\gn> = 0. Notice 
that the function \1/ in (2.10) is defined up to addition of an arbitrary harmonic function 
respecting the boundary condition, and satisfies certain additional relation to be discussed 
below in detail. 

To conclude this section, we single out two special cases of the equation (2.7). The 
first (stationary case) is when the quantities Q and ^ are time - independent. We then 
have from (2.7): Q = Fq(^), where Fq(^) is an arbitrary function, and 

A^ = F (^), (2.12) 

that is, we obtain an equation of elliptic type in dimension (2+0). Most interesting are 
the completely integrable cases of the equation (2.12) with functions Fq(^/) of the form 
±e ± *, icosh 1 !', isinh^, ±sin\l/, — e^ 2 ^ (elliptic version of the Tciceica's equation) 
and some others, which allow to apply the inverse scattering method. A range of such 
cases, with applications to various physical phenomena, is considered in literature (see, 
for instance, [10-12]). The methods used for solving (2.12) include the Hirota method, 
Darboux and Backlund transformations [13-15], the finite-gap integration [14], and the 
inverse scattering transform both for the problem on the whole plane (x, y) [16] and on 
the half-plane {(x, y) : y > 0}[17]. 

Another important nontrivial partial case of the equation (2.7) arises as follows. As- 
sume that f2 = A\l> = Fi(fy) + g with a certain function F 1: and g = g(x,y,t). Assume 
further that either Fiy = 0, or ^ t = 0. We then obtain from (2.7), 

2 it is easy to see that it can be re-written in terms of the stream function alone. 
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9t = ^y9x ~ ^x9y (2.13) 

. This equation, which is given, for instance, in [18], we shall call the simplified vortex 
equation. As well as (2.7), this equation is of significant physical interest, and, impor- 
tantly, is completely integrable. 



3. EXACT SOLUTIONS OF EQUATION (2.13). 



In this section we construct an exact solution of the simplified vortex equation (2.13) 
assuming that g = A\l/ + f(x,y,t), where f(x,y,t) is a known function. 

It is easy to verify by straightforward computation that the equation (2.13) can be 
represented as a compatibility condition for the following system of scalar equations 3 : 

(g-\)(f = 0, Iff = ^y(p x - ^ x (Py, (3.1) 

where ip = ip(x, y, t) is a complex-valued function and A G C is the spectral parameter. It 
is convenient for our purposes to put the first equation of the pair in a somewhat different 
form to obtain, 

9x<Py = g y <Px, <Pt = ^yVx - Vx<Pv (3-2) 

The existence of such representation allows us to construct a family of exact solutions 

of the equation (2.13). To this end, we use the method of the generalized Darboux 
transformation 4 (see, for example, [19]). 
Let 

Q(y, yi ) /■<***> 

<p{x,y,t) = — -, Q(<P,<Pi)= (3-3) 

Mx,y,t) J(xo,y ,t) 

where q(ip,ipi) = (ip x <pi — (p(p\ x )dx + (<f y <fi — (pi y (p)dy is a 1-form, and <fi is a fixed 
solution of (3.2); then the integral Q(<p, fi) does not depend on the contour of integration 
by virtue of the first equality. Let us verify the covarianse of the system (3.2) under the 
transformation (p — > <p, g — > g, ^ — > ^. For the first equation this leads to a dressing 
relation of the form 

9y = 9x^, (3.4) 

Vlx 

which implies that 5 

g(x,y,t) = fo(<fi(x,y,t)), (3.5) 
where /o(vi (•••)) ^ s an arbitrary function. The verification of the second relation gives 

3 In the present form this Lax pair was suggested by A.Ja.Kazakov. 
4 This method is also known as Moutard's transform. 

5 For a first order partial difTerential equation of the form {Y, A}h — 0, with a given A = A(x, y), the 
function Y = Rq(A(x, y)) is a solution for any function Rq(.). 



4 



Qt - Qx^y + Qy^x + Q[-(lny9i)i + ^(ln ^i)* - ^(ln^y = 0. (3.6) 

Solving this linear first order partial differential equation (recall that the function \& is 
known), we find Q explicitly. 

Taking (3.5) into account, we have from (2.11), 

*(x,y,t) = J l D G (x-x\y~y\t)[f ( Vl (x',y',t))-f(x',y',t)}dx'dy', (3.7) 

which solves the initial problem. 

Clearly, the procedure suggested can be repeated and one can construct the iV-dressed 
solutions, but we shall not dwell on that in detail here. 



4. SOLUTION OF THE VORTEX EQUATION. 



1. Linearized equation. 

Let us first consider the linearization of the equation (2.7). Let \l/(x, y, t) = ^i(x, y, t) + 
^o(x,y,t), Q(x,y,t) = Qi(x,y,t) + f2 (x, ?/,£), where Qq = A^o = 0, is a known 
solution of (2.7). We then have 

(A*x) t = (Att^tto, - (Atti^os- (4.1) 

Taking A^i in the form 



A^i ~ a(a, j3, Co) exp{i(ax + (3y + Lot)}, where a(a, /3, Co) — a(—a, —/3, —u>), (4.2) 

and substituting it in (2.7), we obtain the following dispersion relation for the waves 
(assuming that ^/ 0x = p , ty 0y = q , where p , q are constants) 

u = p a - q p. (4.3) 

2. Certain partial solutions. 

Using various approaches, we now construct a set of partial solutions of (2.7). 
a). Let *(x,y,t) = (X(x) + Y(y))T(t). We then have from (2.7): 

a{Xi x + Yiy) = Xi xx Yi — YlyyXl, 

where X\ = X x , Y\ = Y y , a = T t /t 2 = const. Searching for the functions X±, Y\ polynomial 
in x and y respectively, we find, 

V(x, y, t) = [c(x 2 -y 2 ) +d 1 x + d 2 y + h]^, Q = = 0, (4.4) 

where t > 0, and c, d\, d 2 , h are constants. This solution is quadratic in the spatial 
variables and slowly (sub- exponentially) decreases in time. Obviously, it describes a 
vortex-free flow of the fluid. 
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b). Let V = X(x) + Y(y) + T(x, y, t). We then have, 



(T XX -\~ Tyy)t (X XX -\~ T XX -\- Tyy^ X (Yy ~\~ Ty) ( 1 ^ yy ~\~ T XX ~\~ Tyy) y{^X X ~\~ T X ) . 

Letting X = a^x 2 + b\X + Ci, F = a 2 y 2 + b 2 y + c 2 , T = ccos(o;x + /3y + yt + 5), where 
a,, 6j, q, i = 1,2, c, a, /3, 7, 5 are constants, we obtain: a\ = a 2 = 0, 

V(x,y,t) = b±x + b 2 y + d + ccos(ax + fly + 7* + 5), 

(4.5) 

J) = - c (a 2 + /9 2 ) cos(ax + fly + 7* + 5), 

and 7 = a&2 — flbi- These formulae were used, in particular, in [20] to construct solutions 
on the basis of group-theoretical considerations and describe a stream function linear in 
spatial variables on the background of spatial-temporal oscillations. 

In a similar way, one can obtain a solution growing in the absolute value in all variables: 

ty(x, y, t) = b±x + b 2 y + d + ccosh(aa: + fly + yt + 5), 

(4.6) 

Q = c(a 2 + fl 2 ) cosh(a:r + fly + yt + 5), 

where 7 = ab 2 — flb\. 

c). Let us define the complex variables z = x + iy, z = x — iy, and let d = (l/2)(d x — 
id y ), d = (1/2) (c^ + idy). On passing to these variables in (2.7), we obtain 

V ggt = 2i(V ggg V g - y zz - z ^- z ). (4.7) 
Now, setting x — ^z, we have, 

Xzt = 2i(xzzX-XzzX)- (4-8) 

Assume that Xz = 0, that is, x is an analytic function in the domain D. Then it satisfies 
the following nonlinear equation, 

Xzt = -2%XzzX- (4-9) 
In the simplest case when x( z , t) = Z(z)T(t), this gives, 

^ = a 2 , Ei{\nZ) = Ei{\nZ Q ) + < ^e^{z-z Q ), (4.10) 

where Ei (u) = J" (e* /t)dt is the integral exponential function, a, C\ are arbitrary 
constants, and Z$ = Z(z ). Formally, one can write, 

2 °i 

Z( Z )=e Ei 1 [ln^o+| r e^"(z-zo)]_ (4-11) 

Another solution of (4.9) is found if we plug x i n the form x( z ^) = x(0> where 
£ = az + flt + 7, a, fl, 7 being constants. The solution is then given by 
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X(z,t)=£. (4.12) 

In both cases (4.10) - (4.11) and (4.12), using the Green-Cauchy representation (the d - 
problem relation) [21], we obtain for the function \1/ in the case of a bounded domain D, 

*(M,*) = — / Xoi y ,t) dC + — [ I ^iM^AdC, (4-13) 

where the domain D is assumed bounded, and <9D stands for its boundary, x^[z,z,t) = 
V(z,z,t)\ zedl0 . 

2. Lax representation and the method of Darboux transformation. 
The equation (2.7) can be represented as a compatibility condition for the following 
over-determined system of scalar equations [4, 22] 6 : 

L<p = \ip, ipt = -Aip, (4.14) 

where Lp = {Q,p}u, Ap = A G C is the spectral parameter, and the symbol 

{■■}h refers to the "hydrodynamical" bracket defined by (2.6). Here the compatibility, 
and hence the integrability of (2.7), is understood in a "weak" sense, that is, under the 
additional condition 



{Q t + {^,n} H ,p} H = 0. (4.15) 

We now apply the method of Darboux transformation to construct a solution of the 
equation (2.7). Two different versions of the method will be considered. 

i). We use a matrix transform. It is easy to notice that the system (4.14) can be 
rewritten in a matrix form (taking into account that ^ = ^f, Q = Q), 

n x $y - n y <5> X = ^A, $ t =Vy$ x -V x $y, (4.16) 

where 

$ = $(x,y,t,\) = ^° , A = diag(X, A), 

and the condition (4.15) is assumed to hold. 

Let us check the covariance of the equations (4.16), $ — > $, * — > ^, f2 — > fi, under 
the transformation 

$ = $-o"i$A, <7i = ^Ai^ 1 , (4.17) 

where <3>i = $(x,y,t, Ai) is a solution of the system (4.16) corresponding to a fixed value 
A = Ai. We then arrive at the following "dressing" relations, 

= = {Q x ,n} H = 0, {n,a 1 } H = 0, Cl x = Q x , (4.18) 

from which we infer that 



6 From the viewpoint of the spectral theory the first equation in (4.14) can be interpreted as an 
eigenvalue problem for the operator of "generalized" rotation of the plane, Q X J^ —£l y J^, which is easily 
shown (sec, for instance, [2]) to be antisymmetric. 
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* = tf + i?i, ft = ft + Gi, (4.19) 

where iq = Fi(t), G± = Gi(t) are arbitrary functions of time. 

Notice, further, that the substitution * -> * + iq(t) + h(x,y,t), ft -> ft + Gi(f), 

h = h(x,y,t) =^2<U(t)hi(x,y), (4.20) 

i=l 

where the integer L and the functions dj(i), i — 1, L are arbitrary, and hi(x,y) are 
harmonic, does not change the form of the equation (2.7) if the condition 7 

L 

G w = {ft, X^WMtf, (4.21) 

i=i 

is satisfied. 

This allows to write down the new solution of the equation (2.7) as follows, 

f (x, y, f) = y, f) + Fi(i) + h(x, y, t). (4.22) 

ii). We now provide an alternative Darboux transformation, which apparently makes 
it possible to construct a wider class of solutions, as compared to the above. To this end, 
we check the co-invariance of the system (4.14) under the transformation tp — > <p, — > 

ft — > ft, where 

dp = T(x, y, t)((p x - a x tp), (4.23) 

oq = (lny»i) x , </q is a fixed solution of (4.14) corresponding to A = Ai, and T = T(x,y,t) 
is a function to be determined later 8 . Matching the coefficients at (p, ip x and ip xx , we then 
obtain, after a cumbersome calculation, the following over-determined system of dressing 
relations, 

{ft, ft},, = 0, {# - ft} ff = 0, (4.24a) 

AT[(Tl " ft^lf + ai)] ~ { ^ Tai}ff = °' (4246) 
{ft, lnT}^ + n x [&) x - <jA + ft,^ + \& - 1) = 0, (4.24c) 

X X X 

(Tajt + {§, Ta^H - AT[ ^^ (ft^ + ft ra ) - |p] = 0, (4.24d) 
(lnT) t + * x [(hiT) v + (% x - + * y [-(lnT) x + + 

" "X "CP, 



(4.24e) 



3? 3? 3? 3? 



7 This remark can be generalized to the case h(x, y, t) — J2f=i hi(x, y, t), Ah = A/i, = 0, i = 1, L. 
8 In paper [22] the transformation (4.23) was used in the case A = 0, T(x,y,t) — 1/Sl x . 
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Let us analyze the compatibility of this system, limiting ourselves to simpler reduc- 
tions. Obviously there are two essentially different cases: A = and A ^ 0, when the 
coefficients at A should vanish. 

Let A ^ 0. We then infer that V x = V x from (4.24e); V xx = from (4.24d); Cl x /Q x = 1 
from (4.24c). On taking into account (4.24a), we obtain that 

§ = [a (y,t) + c 1 ]x + c 2) 

where Qq(,) is an arbitrary function, and Ci, c 2 are arbitrary constants. From this we find 
that Q = a 0yy x, Cl y = a 0yyy x, Cl x = a 0yy , and that Q y /Q x = (\na 0yy ) y x in view of the first 
equality in (4.2a). Eventually, we arrive at a rather complicated over-determined system 
of equations on the function a , which is likely to be incompatible. 

Now let A = 0, T = 1. On calculating, we then obtain that (4.24) transforms into a 
system of the form 

{n,n} H = o, {Ci,a 1 } H = o, {^-y,n} H = o, {£i,* x } H = o, 

(4.25) 

(^)x = 0, a lt + {^,a 1 } H = 0. 

Substituting 

* = # + £l = n + V 1 , Vi = AE/i, (4.26) 
in (4.25), we get the following system of equations on U\\ 

{n + v u n} H = o, {n + v 1 ,a 1 } H = o, {n + v 1 ,* x + u lx } H = o, 

(4.27) 

{U 1 ,Q} H = 0, <tu + {* + U 1 ,<t 1 } h = 0. 

The compatibility of the first three equations becomes obvious if we set Q + V\ = 
p (Q,<Ji,^/ x + Ui x ), where Po(-) is an arbitrary function. It then follows from the first 
relation that {AUi,Q}h = 0, that is, AU\ = qo(£l) where qo(-) is an arbitrary func- 
tion, hence U\ = A~ 1 qo(Q) = go(ty- This means that the fourth equation in (4.27) 
is compatible with the first three. Moreover, the fourth equation in (4.27) implies that 
U\ y /U\ x = Q y /Q x , and the requirement of the compatibility with the remaining equation 
gives 

Vl ( X ,v,tM=- r °",;" i r li i. (4.28) 

under the assumption that U\ = Ui(x, y, t, Ai), Ui(— oo, y, t, Ai) = and {Q,ai} H ^ 0. 

Thus, we have obtained dressing relations allowing to construct a new explicit solution 
of the equation (2.7) from a known one. It is also possible to obtain the corresponding 
multiple dressing formulae in the spirit of the Darboux transform method [23] . 
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5. CONCLUSION 



Several procedures for constructing the solutions of the 2D vortex equation have been 
suggested in this paper. A number of questions remains open, however: the meaning 
of the notion of " weak integrability" , the possibility of application of the multidimen- 
sional inverse scattering transform method, the Hamiltonian formalism, and, especially, 
quantization of the vortices, to mention a few. 

Notice that various extended versions of this equation are of doubtless interest. For 
instance, the equation for the so-called Rossby waves [6, 7] has the form 

Q t = {Q,^} H -(3Q, 

where (3 is a parameter, that is, it differs from (2.7) by an additional term. As shown 
in [4], it also admits a Lax representation (which is easily recovered from (4.14)), hence 
exact solutions can be constructed. 

Let us emphasize that there is a wide range of adjacent topics related to the theme of 
this paper: the problem of decay of a turbulent spell, models of vortex chains, behaviour 
of vortex ensembles in geophysics, plasma physics, and many others. These reasons urge 
a further analysis of the equation (2.7) and its generalizations. 

The author is grateful to A.V.Urov and P.P.Kulish for support. 
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